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We consider the measurement of the third cumulant of current fluctuations arising from a point 
contact, employing the transitions that they cause in a quantum detector connected to the contact. 
We detail two generic detectors: a quantum two-level system and a harmonic oscillator. In these 
systems, for an arbitrary relation between the voltage driving the point contact and the energy scales 
of the detectors, the results can be expressed in terms of an effective detector temperature T c g. The 
third cumulant can be found from the dependence of T e ff on the sign of the driving voltage. We find 
that proper ordering of the fluctuation operators is relevant in the analysis of the transition rates. 
This is reflected in the effective Fano factor for the third cumulant measured in such setups: it 
depends on the ratio of the voltage and an energy scale describing the circuit where the fluctuations 
are produced. 



I. INTRODUCTION 



The statistics of current fluctuations in mesoscopic 
conductors have been at the center of interest within the 
last decade or so. This statistics can be described by the 
characteristic function, which is the Fourier transform 
of the probability density for a given value of the cur- 
rent. Characteristic functions for many different types of 
systems have been calculated, ranging from the simple 
case of a tunnel junction with Poisson distributed cur- 
rents, to point contacts with binomial distribution and 
to more complicated systems composed for example of 
superconductorsi^ Along with the characteristic func- 
tion, the fluctuations can be described by the cumulants 
or moments of the distribution, the latter being of the 
form 



M n (ii,...,g = ( Yl si(u)), (i) 

i— l,n 



where 8I(ti) = I(ti) — (I), I(t) is the instantaneous value 
of the current, and the brackets (•) refer to ensemble aver- 
aging. Typically described observable is then the Fourier 
transform of M n (Q, t2—ti,...t n — t\) with respect to the 
time differences t» — ij. , and often only the limit where all 
the frequencies are taken to zero is known. 

As many of the studied systems require quantum me- 
chanics for their description, a natural question is the 
proper generalization of Eq. Q to include the fact that 
the current operator I(t) may not commute with itself at 
different times. The answer to this question depends on 
how the fluctuations are to be measured. Levitov, Lee 
and Lesovik suggested the use of a spin coupled to the 
fluctuating current as an imaginary detector £ More pre- 
cisely, when the coupling is turned on at time t = 0, the 
angle of the spin coupled only to the fluctuating current 
I(t), with no average fields, starts to precess along with 
the current. Then the angle at a later time t, averaged 



over the fluctuations, is 
a+(t) = <r+(0) 

x (fexp^i^ £ dtI(t)jTexp^-i^J^ dtl(t)j\ . 

(2) 

Here T and T denote the time- and anti-time-ordering 
operators and g is the coupling constant. This combina- 
tion of time-ordering operators is also called the Keldysh 
ordering. For a classical current, ignoring noncommuta- 
tivity, Eq. J2J would yield the characteristic function of 
the charge Q = f Q I(t) transmitted in a conductor in a 
given time t. Hence, Eq. J5J is one possible generaliza- 
tion of the characteristic function for a quantum current 
operator. In the language of quantum two-level systems 
(TLS's) or qubits, Eq. (J2J describes the dephasing of a 
non-biased (zero level spacing) qubit. 

The advantage of the spin detector is the fact that the 
Keldysh-symmetrized characteristic function of a point 
contact with transmission probability T n for channel n is 
a product of binomial characteristic functions, described 
by the probability T n of success. For example, the first 
three moments of the distribution of currents measured 
in that way are 

(I) =G VY / Tn 

71 

(51 2 ) = G Q eV ^T„(l - T n ) = F 2 e(I) 

n 

(51 3 ) = G e 2 vY,T n {l - T n )(l - 2T n ) = F 3 (I). 

n 

Here Go = e 2 /h and the summation goes over the spin 
and the different transverse channels in the point con- 
tact. Such a characteristic function thus has a transpar- 
ent classical interpretation. 

But spin precession or qubit dephasing in real time 
is difficult to measure. Also, Levitov-Lee-Lesovik spin 
detector responds only to the zero- frequency noise. An 
alternative method is to measure the excitation and re- 
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laxation rates of a quantum system subject to a fluctuat- 
ing force produced by the current. This is the approach 
taken in the present paper. These transition rates deter- 
mine the static state of the density matrix for the system, 
which can in many cases be described by an effective tem- 
perature. Measuring the effective temperature is much 
more simple than spin precession, and it also gives ac- 
cess to the finite-frequency moments of the noise. In this 
paper we discuss how two simple quantum detectors be- 
have in response to a non-Gaussian fluctuating force. An 
exact solution of this problem with the full characteristic 
function of fluctuations is not known to us, and there- 
fore we resort to expanding in the cumulants of these 
fluctuations up to the third cumulant. 

It turns out that in this case the proper way to define 
the third-order correlation function is by time-ordering 
two of the current fluctuation operators, and leaving the 
third one free (c.f., Eq. (J3J)). 4 The relaxation of the qubit 
in this case is sensitive to the frequency dependence of 
the third cumulant^ and the relevant Fano factor "mea- 
sured" in this process depends on the relation between 
different frequency scales of the problem: If the qubit is 
"close" to the noise source, so that the frequency depen- 
dence is given by the voltage V applied over the scat- 
terer, the third cumulant strength is characterized by 
F 2 - F 3 = 2 £„ T„ 2 (l - T n )/ T n . This factor is very 
small for a tunnel junction, but finite for other types of 
scatterers. However, placing the qubit further from the 
noise source introduces another frequency scale into the 
problem, characterizing the circuit between the source 
and the qubit. For voltages above this frequency scale, 
the third cumulant effect on relaxation is characterized 
by the "usual" Fano factor, F%, In this limit one repro- 
duces our results in Ref. 0, where we assumed that the 
frequency dependence is solely governed by the circuit. 
Besides the Fano factor, also the dependence on the volt- 
age V is different in the two limits: in the first case it is 
logarithmic and in the second case linear. 

Alternative "on-chip" detector schemes for measur- 
ing non-G aussian fluctuations have been suggested in 
Refs. IgOt MoI 1(11 ill 12l and a few schemes have already 
been experimentally realized, see Refs. IT^Il4lll5lll6i 



sufficiently weak, i.e., KT env *C f2, such that it can be 
taken into account perturbatively. The detector response 
will depend on the ratios between the voltage and the 
detector level spacing, eV/Q (or, at a finite temperature, 
on fcsT/rj) and on the ratio Huj c /Q. Finally, the response 
depends on the relative magnitude of the transition rates 
coming from the noise source and of those coming from 
the Gaussian bath. 
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FIG. 1: Schematic idea of the fluctuation measurements: Cur- 
rent fluctuations in the noise source are coupled to a quantum 
detector, where they induce transitions between the detector 
energy levels. Another source for the transition rates is the 
intrinsic Gaussian environment of the detector, which can be 
modeled via an ensemble of harmonic oscillators. Up to the 
third order in the coupling coefficient g, the transition rates in 
the detector depend on the second and third cumulant of the 
fluctuations. As the latter is odd in the driving average cur- 
rent through the noise source, its effect can be measured by 
detecting the change in the rates when the sign of the driving 
average current is reversed. 



II. EFFECT OF THIRD CUMULANT ON AN 
ARBITRARY DETECTOR 

Consider the system depicted in Fig. ^ A non- 
Gaussian noise source is coupled to a detector whose state 
we aim to describe. The Hamiltonian of the system can 
be decoupled into 



A. Frequency scales and different regimes 

The noise source - detector system can be character- 
ized with a few frequency/energy scales whose relative 
magnitudes determine the detector response. The main 
frequency scales are that coming from the detector level 
spacing Cl/ti, one characterizing the circuit connecting 
the noise source and the detector, say oj c (for an exam- 
ple, see Ref. 0), and the rate r cn v of transitions in the 
detector induced by its own environment. Furthermore, 
the fluctuations in a point contact are characterized by 
two further scales, given by the temperature T and the 
voltage V over the contact. 

In this work, we assume that the level broadening is 



Ft — -f^dct "f" H env -\- Hecoup ~t~ F[ no [ sc -f" H COU p . 

Here -ffdet is the Hamiltonian of the detector, specified in 
Sees. II VI and IVl It is in general described through a set of 
collective variables with mutually noncommuting opera- 
tors Ai . This detector is coupled to its own environment 
described by H cnv through the coupling H ccoup . We as- 
sume this environment to be Gaussian, such that we can 
model it via an ensemble of harmonic oscillators^ 

3 

Recoup = 1£ ^ ^ + ^ + f^ A2 12 A '' 
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where a(t> i S the bosonic annihilation (creation) operator 
of an oscillator in the bath, u> n its eigcnfrcqucncy, and Xj 
is the coupling constant from this oscillator to the system 
variable A. 

Finally, the noise source is described by i? no isc and 
coupled to the detector via H coup . The latter connects 
in general a set of collective variables of the noise source 
to another set of variables of the detector. As we aim to 
describe the measurement of current fluctuations, we ex- 
plicitly assume that the previous is the current operator 
I in the noise source^ Thus the coupling is of the form 

e t-r 1 

i 

where gi are dimensionless coupling constants. We fur- 
ther assume that I commutes with A4. 

Including the effects of H coup up to the third order in 
the coupling gi, we do not have to specify H nomc , but 
it suffices to concentrate on the different correlators (cu- 
mulants) of I. For simplicity, we include the effect of the 
average current on the detector Hamiltonian iJdot, such 
that we can take (I) = 0. Here the brackets (•) denote 
quantum averaging over the density matrix of the noise 
source. In particular, we concentrate on the noise power 
spectral density, 



dte luj{t - to) (i(t)I(t )), 



and the partially time-ordered third cumulant^iS 

/>oo 

8 3 I(uj u u 2 ) = / d(h - t )d(t 2 - t )e lWl(tl - to) 



x e^^-*°)(f [7(t )/(ti)]J(t 2 )). 



(3) 



In a static system considered in this paper, these corre- 
lators are independent of the time to- 

Assume the uncoupled detector is described via the 
energy eigenstates |n) with energies E n , i.e., H^ e t\n) = 
E n \n). Without coupling to the environment, the state 
of the detector is described by a density matrix 
which in general may show coherent oscillations between 
the different states. Assume now we turn on the cou- 
plings i/coup and -ffecoup at some time t . If H coup 
or -Hecoup do not commute with -ffdet, after some time 
Pnm(t) tends into a diagonal steady state formA These 
diagonal entries P n = pjf£(t 3> to) are obtained from a 
detailed-balance relation of the form 



P F 



(4) 



In addition, the total probability has to be conserved, 
i.e., ^2 n P n = 1- Here r m _>„ is the total transition rate 
from the energy eigenstate m to the eigenstate n, due to 
the coupling to the environment. 

Up to the third order in the coupling constants gi and 
Aj, the transition rates originating from the coupling to 



the oscillator bath and the noise source are uncorrelated 
and we may write 



pcnv 1 pnoisc 

>ri m — >n ' m — >n' 



Here 
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£*3 



JlUJmn / (fcT e )penV 



is the transition rate from state m to state n due to 
the coupling of the detector to its Gaussian bath with 
temperature T e , up to the second order in the cou- 
pling constants Xj. Here A mn = (m\A\n) and io mn = 
[E m — E n )/h. In what follows, we use a short-hand no- 
tation A 2 = J^j ^j- The transition rates induced by the 
noise source are of the form^S 

rr?„ = i^ 5l 2 iAri 2 5z(^ mn )+r^L„. 



The rate from the third-order term is 



r (3) 



du> 



-UJ„ 



LU - LU hl - IT) 



A m ^ A ^ n A nm 
" /I,- 



(5) 

Here ry is a positive infinitesimal. In this paper, this 
integral is evaluated for two generic detectors in the case 
of noise originating from a point contact. 



III. FREQUENCY DEPENDENT SECOND AND 
THIRD CUMULANTS OF A POINT CONTACT 

A general scheme for calculating arbitrarily ordered 
frequency-dependent current correlators from the scat- 
tering theory was laid out by Salo, Hekking and Pekola 
(SHP) in Ref. [3. Their results are applied here in order 
to calculate the response of our generic detectors to the 
third-order current fluctuations. SHP decompose the op- 
erator describing the current through a given scatterer to 
"in" and "out" parts, I = I- m — I ou t, and show that time 
ordering between two current operators, I{t\) and Ifo), 
can be expressed in terms of ordering between L m and 
J out . For the latter, they find that a time-ordered prod- 
uct of a pair of (in, in) or (out, out) operators is the same 
as the unordered pair, and that time ordering a pair of 
(in, out) operators corresponds to an ordering where 7 ou t 
is placed to the left of I ia . 

A practically important noise source is a point con- 
tact with an energy independent scattering matrix. A 
"point contact" in this case refers to a system through 
which the electron time-of-flight tv is much smaller than 
the other time scales of the problem. In this case, 
the second-order correlator can be written in the form 
Si(uj) — Sf(ui) + Sj xc (uj), where the vacuum fluctua- 
tions are Sf(u>) = 2hcuG9(ui), and the excess noise is 
given by^ 1 
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fhuj\ \ eysinh(f£)-2ftwcoth(^)sinh 2 (^) 

<Sf (w) = Ghw coth ( — — J - sgn(w) ] + F 2 G- 



V2fcTy 6 v 7 cosh(f^) - cosh(ff) 

F 2 G(e\V\ - h\w\)0(e\V\ - h\w\). 



Here G is the conductance of the point contact, V is the 
voltage applied over it, and F 2 = E n T n (l-T„)]/ ^„ r n 
is the Fano factor characterizing the transmission eigen- 
values T n of the contact. Written in this way, the excess 
noise is a symmetric function of frequency, and thus it 
contributes to excitation as much as to relaxation. 

For the partially time-ordered third cumulant in the 
case of a point contact, one can deduce from the results 
of SHP 22 



<5 3 /(wi,w 2 ) = S ioo (-Ui - lu 2i -lu 2 ) + S ioo (u>i,-w 2 ) 

+ S ioo (-UJ 2 , -U>x - LU 2 ) - S 000 (—Ul - W2 5 — W2). 

Here 



(6) 



S ioo (u;i,uj 2 ) = eF 2 G[A(uJi,u} 2 - v) - B(ui,wi - uj 2 - v)] 
Sooo = eF 3 G[A(uji , lu 2 - v) + A{w\ - v, lu 2 ) 

+ A{uj\ + V, LU 2 + v) — B(jjL>±, 0J\ — lu 2 — v) 

— B{ii)\ — V,UJ\ — UJ 2 ) — B{uji + V, W\ — u 2 + v)], 

v = eV/Ti and F 3 = ^ n T n (l - T„)(l - 2T n )/£ n T n . 
The functions A(uji) and B(lj 2 ) are defined as 



A(x u x 2 ) = J dEf(E)(l-f(E + Xl ))(l-f(E + x 2 )) 
9(xi)9(x 2 ) mm(xi, x 2 ) 

(7a) 



B{x u x 2 ) = J dEf(E)(l- f(E + Xl ))f(E + x 2 ) 



(7b) 



(xi)8(x 1 — x 2 ) min(xi, X\ — x 2 ), 



where f(E) is a Fermi function. 

In an electric circuit containing reactive elements, the 
fluctuation spectra are modified in a frequency dependent 
way. These types of modifications can be fairly generally 
calculated with a Langevin approach (see for example 
Refs. |2j| and |3 f° r an exception relevant for the third 
cumulant, see Ref . l25fl . In the case of a time-independent 
average current, the noise spectra are modified according 
to 

Sf(u;) = G{lj)G(-lj)SF c (uj), 
5 3 I c (u>i,uj 2 ) = G{ui)G{u) 2 )G{-u) X - u 2 )5*I{u) U w 2 ), 

(8) 

where G(u>) is a function characterizing the circuit. In 
what follows, we choose 



G(u) = 1/(1- iw/w,-), 



(9) 



typical for a circuit with reactive elements next to the 
point contact. 

Below, we aim to calculate the outcome of these spec- 
tra on two generic detectors coupled to the noise source, 
a quantum two-level system and a harmonic oscillator. 



IV. QUANTUM TWO-LEVEL SYSTEM 

The general Hamiltonian for a quantum two-level sys- 
tem is 



H 



B, 



B , 



TLS 



where <r z / x are Pauli matrices, and B z / X the effective 
magnetic fields. We assume that such a system is coupled 
to the noise source via a coupling Hamiltonian 



#cou P = -g{I + SI)a z . 

The field pointing to the y-direction would not add any 
more generality to our model. The average current I 
can be included in the classical control field by defining 
B z = B z + 2hgl/e, so we can only concentrate on the 
fluctuations 51. 

The ground and excited states of this system are given 
by |0) = -I3\ T) + a\ |> and |1) = a\ |) + 0\ I), with a = 
cos(0/2), (3 — sin(0/2) and 4> — arctan(£? z / B x ). The en- 
ergies of these states are _E /i = CI = y/B%+ B%. 

Now, the second-order contribution to the excitation 
rates is 



(l I - ^2 |cr° 1 | 2 S'noiso(-^/?i), 



^(2) 



(10) 



and the third-order contribution can be obtained from 



.(3) 
- 0->l 



4K ol | 2 cos(0)Re 



6 3 i{uj,n/h) 



duj- 

LU + il/n - IT] 



du- 

00 lu-it] 



(11) 



Here we used the fact that ct] 1 = — a z ° — cos(0). Note 
that the matrix element IcrJ 1 1 2 — sin 2 (</>) is common for 

(2/3) 

both rates. The corresponding relaxation rates T\_!^ can 
be obtained from the excitation rates with the substitu- 
tion CI — > —CI. 
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In the case of a point contact at a vanishing temper- 
ature (kgT <C eV,ft), the second-order excitation and 
relaxation rates are given by 



^(2) 



F (2) _ 
1 l->0 — 



2 G 



(e|Vj -n)0(e\V\ - 0) 

1+ (fi/^c) 2 

F 2 (e|y| -n)0(e|V| - 0) 
1 + (f7/?iw c ) 2 



(12a) 
2fi" 

(12b) 



The calculation of the integrals required for the third- 
order effect in this case is detailed in Appendix ^] The 
general result with an arbitrary ratio between eV and 
hu) c can be found analytically, but it is too long to be 
written down here. In the limit Q < \eV\ -C huj c we 
obtain 



.(3) 
0->l 



9 i 01 I 2 



^9 I 01 1 2 



cos{<fi)Gsgn(V){F 2 - F 3 ) 



fill! 



|eV| - fi 



leF I In 



\eV\ 



\eV\ — ft , 



(13) 

and = — Tq 3 ^. For \eV\ < f2, the contributions to the excitation rates from the third cumulant vanish, but 

there is a contribution to the relaxation rate, 



lf_i (|eV| < O) 



9 \affcos{cp)GF 2 



8nh 3 uj 3 n 



4nV + 5nV 2 fi 2 + 4 



V. 



(14) 



This contribution becomes small in the limit Tilj c fl. 
For \eV\ > fl, these rates tend to ±r^ se with 



r^L = ^k z °Tcos(0)Gsgn(VO(F 2 - F 3 )fi 



In 



r 



\eV\ 

n 



9. 



2\eV\ 




(15) 



Note that for a tunnel junction, F 2 = F3, and there is no 
contribution from the third cumulant to the transition 
rates. For other types of contacts, the rates are deter- 
mined according to F 2 - F 3 = 2 T 2 (l - T n )/ T n . 
This is the same as the zero-frequency Fano factor one 
would get for an unordered third cumulant^ 

In the opposite limit where the voltage by far exceeds 
the scale huj c /e set by the circuit, the frequency depen- 
dence is governed by uj c and the third-cumulant effect on 
the excitation rate is given by 



r (3) _ _ 



4irg 3 



,01|2 



cos((/>)F 3 



Gh 3 ^ 



5/iVn 2 



-V 



,(3) 



ft 4 

(16) 

and again = — Tq^. This is the result one would 

obtain by assuming S 3 I c (uji 1 uj 2 ) ~ G(uji)G(uj 2 )G(— lo\ — 
LO2)S 3 I(0, 0) as was done in Ref. 0. 

The third-order contributions to the relaxation rates 
induced by the noise source are plotted in Fig. [21 for a 
few example cases. 

Now let us analyze the detection of noise with the two- 
level system via the steady-state occupation probabilities 
Pq and Pi = 1 — Pq of the states |0) and |1). These satisfy 
the detailed-balance condition, Eq. (0J . Analogous to the 
equilibrium system, we can define an effective tempera- 



ture of the quantum two-level system, 



hi 



(ft) 



hi 



r (2) , r <3) 

r (2) 7 r (3) 
L 0— »l" t_A 0— *1 



+ o(g 4 ). (17) 



The contributions from the second and third cumu- 
lant can be separated by considering what happens to 
Ttls when the voltage across the point contact is re- 
versed. We define the average temperature Ttls = 
(Ttls(V) + 2tls(-V))/2 and the difference AT TL s = 
Ttls — Ttls(— V). In the lowest order in the coupling 
constant g, the previous is then independent of the third 
cumulant, whereas the latter is directly proportional to 
it. 

The limiting case expressions for Ttls aud ATtls 
depend on the relative strengths of the relax- 
ation/excitation rates from the bath and from the noise 
source, and on the magnitude of the circuit frequency 
scale uj c . The previous is easiest to characterize through 
the ratio of the differences between relaxation and excita- 
tion they cause (this is essentially the " friction" strength 
in classical models)^ 



A = A 



i (^g + » 2 )sin 2 ((/.) 
2h 2 g 2 Guj^n ' 



(18) 



For fl <C ksT e , \eV\, we then get an asymptotic expres- 



FIG. 2: (Color online): Third-order contributions to the ex- 
citation (solid lines) and relaxation (dashed lines) rates of a 
quantum two- level system coupled to a point contact. The 
rates are plotted for four different types of contacts with 
equal conductance G: tunnel contact (blue, F2 = F3 = 1), 
dirty interface^ (black, F 2 = 1/2, F 3 = 1/4), diffusive wireS 
(red, F2 = 1/3, F3 = 1/15), and a chaotic cavity— (green, 
F2 = 1/4, i<3 = 0). The other parameters used in this plot 
are <f) — 7r /4 and T\w c — 20Q. Inset shows the low- voltage 
region. The rates are proportional to the dimensionless con- 
stant GRk where Rk = h/e 2 is the resistance quantum. 



sion for Xtls, 

2k B T e A + F 2 \eV\ (1 - F 2 )fi 



/ 



2(1 + A) 



2(1 + A) \k B T e 



(19) 

For |eV| 3> k B T e ,fl, there is another asymptotic expres- 
sion, 



k B F 



B^TLS 



F 2 \eV\ 1- J 1 2+Acoth(^^j ^ / q 
2(1 + A) + 2(1+ A) \\eV\ 

The latter equation is valid for an arbitrary ratio between 
f2 and kT e . 

The effective temperature Ttls as a function of the 
voltage V is plotted in Fig.[31for a few types of junctions 
and in Fig. 0]for a few values of uj c . 

The antisymmetric part ATtls of the temperature 
with respect to the voltage through the scatterer depends 
strongly on whether the frequency dependence is gov- 
erned by the voltage or by the circuit. In the previous 
case, for fi <C k B T e , \eV\ <C Tilu c (noise source "close" to 
the detector), the asymptotic expression for ATtls is 



k B ATtls = sgn(]/) 



2.g(T 2 - F 3 ) cos(^) 



x (F 2 \eV\ +2k B T e A) 



In 



(1+A) 2 
eV\ 



n 



n 



(21) 



For a noise source placed "far" from the detector, i.e., 



FIG. 3: Effective temperature of the two-level system as a 
function of the voltage through different types of point con- 
tacts coupled to it. The different curves correspond, from top 
to bottom, to a tunnel contact (blue), dirty interface (black), 
diffusive wire (green) and a chaotic cavity (red). The behavior 
is mostly dictated by the Fano factor T2. Other parameters in 
the plot are g = 0.05, hco c = 20Q, <f> = tt/4, A = 1, k B T e = fi, 
and G = 1/R K - 





FIG. 4: Effective temperature of the two-level system as a 
function of the voltage through a diffusive point contact cou- 
pled to it via circuits characterized with different u c . In the 
right, from top to bottom: hto c /Q = 100 (green), 10 (red), 
0.1 (blue) and 1 (black). Other parameters are as in Fig. [2] 



f2, hu c <§; \eV\, kT e , we get 



k B ATtls = -^cos(<£) 
x (F 2 \eV\ + 2k B T e A) + o 



AgF^TThuJc 



Q, TiLj c 



■ Ay 



(22) 



These limits arc illustrated in Figs. [S] and which show 
the temperature difference ATtls as a function of the 
voltage over the point contact for different types of point 
contacts and for different uj c , respectively. For a tunnel 
junction, ATtls/^ is negative for all values of the volt- 
age as the logarithmic term in V (Eq. I|21|l ) is absent, 
whereas for a chaotic cavity F3 — and the absence 
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of the linear term (Eq. I |22| l) leads to a positive defi- 
nite AT TL s/V. For other types of junctions, there is a 
crossover from positive ATtls /V to a negative ATtls /V 
as eV roughly crosses %u> c . 




V. HARMONIC OSCILLATOR 

Below, we discuss two schemes for using a quantum 
harmonic oscillator as a detector of the current fluctu- 
ations. The first scheme couples a displaced oscillator 
to the fluctuations through the position operator x , and 
the second scheme through the second power of the mo- 
mentum operator, p 2 . It turns out that at least up to 
the third order in the coupling operators Aj, the first 
scheme is insensitive to the third cumulant, whereas the 
latter follows quite closely the qubit scheme. First, we 
note that a simple harmonic oscillator coupled to the 
fluctuations via its position operator, i.e., A = x directly 
implies r' 3 ) = as Eq. J5J requires that all the matrix 
elements between the initial, the final, and one or more 
intermediate states are nonzero. This cannot be satisfied 
as x couples only the neighboring states. 



FIG. 5: Asymmetric part of the temperature ATtls as a func- 
tion of the voltage over the point contact for different types 
of point contacts. In the right, from top to bottom: chaotic 
cavity (red) , diffusive wire (green) , dirty interface (black) and 
a tunnel junction (blue) . Other parameters are as in Fig. [3] 
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FIG. 6: Asymmetric part of the temperature ATtls as a func- 
tion of the voltage over a diffusive point contact with different 
frequencies characterizing the circuit: in the right, from top 
to bottom: huj c /tt = 100 (green), 10 (red), 0.1 (blue), and 
I (black). Other parameters are as in Fig. [3] The curve for 
fibi c — 100S1 is in the logarithmic regime, Eq. 12H . for all 
plotted voltages, and it is characterized by the Fano factor 
F-2 — F3. The curves for hui c < 1 are in the linear regime, 
Eq. 12211 . and the case hu c = 10fi has a crossover between 
the two. Inset shows the low-voltage regime, reflecting the 
behavior described in Eq. I|f4|l . 



The relative temperature change ATtls/^tls is thus 
logarithmic in V for n,ksT e <C \eV\ <C Tilo c and linear 
for \eV\ > Q,,k B T e ,huj c . 



A. Displaced harmonic oscillator coupled to the 
position operator 

This restriction may be overcome by considering a dis- 
placed harmonic oscillator, so that the effective coupling 
is to x — Xq, where xq is a scalar displacement. However, 
one can quite generally show that even in this case the 
resulting T^ 3 ) = 0. 

To be specific, consider a harmonic oscillator coupled 
to the fluctuating current. The Hamiltonian for the os- 
cillator is 



h 

2m 



Assume the fluctuating current is coupled to the in- 
coordinate, i.e., the coupling is described by 



-(lb + Si)\/2mhujQx' . 



Due to the average current term lb, the oscillator po- 
tential minimum is shifted from x — to xq = 
2hglb l(&\J ttlluq ) . Defining a new displaced operator 
- xq we get, neglecting the unimportant scalar 



x — x 
terms, 



Hi 



HO 



H 



2a2 



coup 



2m 



-mu> n x 



+ g 01 (x + x ). 



(23) 



In what follows, we assume that lb can be controlled sep- 
arately from the current flowing through the noise source. 
This type of a separation of average and noise currents 
was discussed for example in Ref . |3(J 

We proceed in the usual way by defining the harmonic 
oscillator annihilation and creation operators, 

I A 1 

(mujQX T hdst) ■ 



\J2hm,LOQ 
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Now the energy eigenstates |n) of the "average" Hamilto- 
nian are those of the number operator, N\n) — a)a\n) = 
\n). The fluctuations are coupled to the operator A = 



\/2mTiuj${x-\-x$). This has a finite matrix element A n>n ±i 
between the neighboring states due to the operator x. 
Moreover, the displacement xq makes the diagonal ma- 
trix element A nn also finite, and independent of the level 
index n. Similar to the two-level system, we can write the 
second-order excitation and relaxation rates due to the 
external noise from a point contact at ksT <C | eVj , hojQ, 



to the transition rates due to the current fluctuations at 
T = are 



.(2) 
n — >n+2 

r ,2 



= i-( n + l)(n + 2)Si(-2u l0 ) 



r(n + l)(n + 2)GF 2 {\eV\ - 2hoj Q )0(\eV\ - 2huj Q ) 



(27a) 



r ( ,^ n = ^(n+l)(n + 2)S I (2u J o) 



■(n+l)5/(-wo) 



(24a) 



= ^(n + l)GF 2 (\eV\ - hw o )0(\eV\ - tiu ) 

r^„ = ^(n + i)S/M 

= ^(n + l)G[F 2 {\eV\ - htu )8(\eV\ - htu ) + 2hoj }. 

(24b) 

The third-order contribution to the excitation rate is 



.(3) 
n — ►n+l 



^(n + l)Re 



du> 



5 3 i(uj,n/h) 

u) + Q/h — ir\ 

6 3 i(u,n/h) 

uj — irj 



(25) 



and the relaxation rate T 



(3) 

n+1 — >n 



can be obtained by re- 



placing luo by — u>o inside the integrals (but not in the 
prefactor). 

These integrals are the same as for the quantum two- 
level system, Eq. Ijllll . up to a sign between them. But 
as shown in Appendix [S] the two integrals give exactly 
the opposite contribution under quite general conditions, 
and therefore I^ 3 ) vanishes. 



= ^L(n + l)(n + 2)G[F 2 (\eV\ - 2huj )6(\eV\ - 2hto ) + 4fwj ]. 

(27b) 

The third-order contribution to the excitation rate is 

„3 



.(3) 
■ n— m+2 



9" 



(n + l)(n + 2)Re 



(2n + l) 



6 3 I(lu, 2uj /Ti) 
lu + 2u>o — irj 



+ (2n + 3) / duj 



5 3 I(uj, 2uj /h) 



uj — irj 



= - : ^r(n + l)(n + 2) [(2n + l)/ 3 + (2n + 3)/i] 

(3) 

and the relaxation rate T n ^_ 2 ^ n can be obtained by re- 
placing wo by —luq. Using the fact that ^3 = — 1\ (see 
Appendix [5), we get 



r 



(3) 

n — m+2 



■^-(n+l)(n + 2)Ji(2w ). (28) 



The third-order rate has thus the same level-dependent 
prefactor as the second-order rate. Therefore, we can 
again define an effective temperature, which now is of 
the form 



2hojQ 



In 



r (2) T r (3) 



(29) 



B. Coupling to the square of the momentum or 
position operator 

Assume one could vary the mass of the harmonic os- 
cillator via the fluctuating current. In this case, the cou- 
pling to the fluctuations would be of the form^i 

iJ coup = S -^lf = -gSI(tf - a) 2 . (26) 



2m 2 



This is independent of the level index n, provided the 
fluctuations coupling linearly to x or p can be neglected. 

Because of the similar form of the rate expressions as 
in the qubit case, the behavior of the effective tempera- 
ture is similar to the qubit, provided the prefactor cos(</>) 
in Eqs. (|15I16I21I22I) is replaced by —1 and the matrix 
element a® 1 by (n + l)(n + 2). In the case of coupling to 
x 2 instead of p 2 (i.e., varying the spring constant rather 
than the mass), the only difference is the inverted sign of 
the third-cumulant contributions to the rates. 



Such an operator A = {a) —a) 2 has finite matrix elements 
between next-nearest neighbor energy levels n and n + 2 
of the oscillator, A n: . n+2 — y/ (n + l)(n + 2) and diagonal 
matrix elements, A n>n = —2n — 1. To be able to further 
describe the system with an effective temperature, we 
assume that the coupling to the bath is much weaker than 
the coupling to the noise source, and the previous can 
hence be neglected. Now the second-order contribution 



VI. EXPERIMENTAL DETECTOR 
REALIZATIONS 

To exemplify the measurement of the third cumu- 
lant through the polarization of a qubit coupled to the 
fluctuating current, we consider three specific exam- 
ples: a persistent current (or a flux) qubit? 3 ^ 3 ^ a phase 
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qubit, 34,35 and a charge qubit. 36 The fluctuation mea- 
surement schemes with these qubits are illustrated in 
Fig. Apart from the measurement schemes, we aim to 
discuss typical values for the coupling constant g, level 
splitting Q, circuit frequency scale u c , and the dimen- 
sionless constant A characterizing the ratio between the 
intrinsic and induced noise. 




FIG. 7: Fluctuation measurement schemes using supercon- 
ducting qubits: a) Persistent current qubit, b) phase qubit, 
and c) charge qubit. 

In all cases, the level splitting fl/h turns out to be of 
the order of a few GHz at minimum. Also, for all systems 
the relevant reactive element causing dispersion of noise 
is the inductance L of the lines feeding the noise current, 
and therefore the frequency scale uj c — \/(GL), where 
G is the conductance of the shot noise source. As L is 
typically of the order of some 10 pH, and R may vary 
between some 10 Q to some 100 kf2, we have oj c ranging 
between 10 3 to 10 6 GHz. This means that unless special 
care is taken to make a large inductance, hio c j£l is at 
least a few hundred. 

However, note that at least in diffusive wires and 
chaotic cavities, there is an additional frequency scale 
given by the inverse dwell time 1/td or the inverse screen- 
ing time 1 /t sc ,2L22iiS For example, for a diffusive wire of 
length L = 1/im and diffusion constant D — 100 cm 2 /s, 



we have 1/td w 10 GHz. When 1/td or l/r sc is less than 
min(|eV|/?i, lu c ), these have to be taken into account sep- 
arately. 

With the qubits, the detection takes place by applying 
a steady current through the noise source, and measur- 
ing the occupation number pn of the higher qubit state 
many times, i.e., averaging over the realizations of cur- 
rent fluctuations. The effective temperature is then 



Pn 



Pn 



The third-cumulant effect can be controlled by tuning the 
phase 4> through the average fields B x and B z . When 
either of these fields vanishes, also the third-cumulant 
effect should vanish. 

The temperature measurement of the harmonic oscil- 
lator depends on its realization: For a true oscillator 
based on a resonant LC-circuit, the temperature should 
be measured either by measuring the current noise power 
in the oscillator, or coupling it to a nonlinear system, say 
a SQUID, and measuring its response. If the realization 
is a current-biased Josephson junction, the temperature 
detection can be done via the measurement of the ther- 
mal escape rate as in Ref. |23 

Voltage fluctuations could also be coupled to the qubits 
or oscillators. Typically these would couple to the per- 
pendicular external field component compared to the cur- 
rent fluctuations. However, the frequency dependence of 
the third cumulant of voltage fluctuations is not known, 
and the response of the qubits might in this case be some- 
what different. 



A. Persistent current qubit 

In the persistent current qubit, a fluctuating current is 
easiest to couple to the qubit current through the mutual 
inductance M as in Fig.[7|i. 32 In the basis defined by the 
clockwise and anticlockwise current, this corresponds to 
coupling to a z . In this case, the coupling constant g is of 
the order of 



9pcq = -MAI qb 



-MI C , 



where A.I q b is the different between the currents corre- 
sponding to the two qubit states. This is of the order 
of the critical current I c of the Josephson junctions. Us- 
ing M « 2 pH and I c ~ 3 /iA close to the experimental 
values we get g s» 0.01. Finally, with the intrinsic 
relaxation time T\ — 1 ms, il/h =10 GHz, R = Rk and 
<f> = 0, we get A as 0.003, i.e., the intrinsic bath effect is 
almost negligible. 



B. Phase qubit 

With the scheme depicted in Fig.[7j3, the external cur- 
rent fluctuations can be again coupled to the diagonal 
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element of the qubit Hamiltonian, i.e., a z 
strength isM 



The coupling 



a) 



9pq 



edEio 
h dl b ' 



where E w is the level separation between the qubit 
states. For the bias current close to I c , this is given by 



huj p (i) 



5 hu) p (I) \ 
~ 36AC/(7) J 



where uj p (I) = 2 1 / 4 ^/8E J E c (l - I/Ic) 1 ^ and AU(I) = 
AV2/3Ej{l-I/I c ) 3 / 2 , Ej = HI c /(2e) and E c = e 2 /(2C) 
are the Josephson and charging energies of the junction, 
and I c is its critical current. Thus we have 



5e c + 3 • 2 3 / 4 Vg^(l 



)5/4 



6(1 -i b y 

where if, = I/I c , and e c = Ec/Ej. With the values^ 
C = 6 pF and I c = 21 //A and i& « 0.99, we get g pq « 0.1. 
In such qubits, the relaxation due to the intrinsic bath is 
slower than w p /1000, which implies A < 0.1. Therefore, 
the intrinsic bath should also here be negligible. 



C. Charge qubit 

Coupling a charge qubit in a form of the Cooper pair 
bo3i2& to current fluctuations takes place via the oper- 
ator <j x in the natural charge basis of the qubit. This 
is accomplished by coupling the fluctuations to the flux 
controlling the Josephson energy of the Cooper pair box. 
Due to this slightly different type of a coupling, the phase 
4> should be defined as = arctan(.R r /.B 2 ), but otherwise 
the rate expressions stay the same. In this case, the cou- 
pling strength is given by^ 

■ 9 = ^kr< 



where M is the mutual inductance and I c (<fr x ) is the crit- 
ical current of the box at the average external flux Q x . 
With M = 2 pH and I c = 300 nA, g can be made to 
vary between almost zero (the "sweet spot" where the 
first derivative of I c (&x) vs. $ x vanishes) to some 10~ 3 . 
A higher coupling strength can be obtained by increasing 
I c and thereby going away from the strict charge basis as 
in Ref. EH 



D. Cooper pair box in a resonant circuit 

Coupling the fluctuations to a flux controlling the 
Josephson coupling of a Cooper pair box placed in a res- 
onant circuit, one may control the effective capacitance 
of the resonator Such a setup corresponds to that 
studied in Subs. IV Bl The coupling strength g depends 




5I(t) 




M 



FIG. 8: Fluctuation measurements with oscillators whose 
mass or spring constant is driven with the fluctuations: a) 
Cooper-pair box and b) driven Josephson junction. In both 
cases, the Josephson energy is tuned with the fluctuations. 
In a), the difference to a charge qubit is in the fact that the 
transistor is placed as a part of a resonant circuit. 



on the relative ratio between the Josephson and charg- 
ing energies, Ej and Ec of the Cooper pair box. For 
Ej = 10Ec, the relative capacitance modulation AC/C 
can be of the order of 0.05 for a change A$ x = h/2e in 
the external flux4 3 » The coupling constant in this case is 



.9 = 



AC e 2 Muj 



For M — 100 pH and ujq =100 GHz, we would hence 
obtain g ~ 10~ 4 . This is quite a small value, and using 
this scheme would require optimization of the Ej/Ec- 
ratio and maximizing both M and uiq. 



E. Current biased SQUID 

Instead of directly controlling the bias current of the 
Josephson junction as in the phase qubit scheme above, 
one may also envisage coupling the fluctuations induc- 
tively to control the flux in a SQUID as in Fig. The 
Hamiltonian of a biased symmetric SQUID is 

Hjj = -AE c dl - Ej cos (^^) cos(^) - ^I b <p, 

where ip is the phase across the SQUID, is the bias 
current, <fr x is the external flux through the loop, and 
$o = h/(2e) is the flux quantum. We assume the self- 
inductance of the loop small enough, so that it can be 
neglected. For If, much lower than the critical current 
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I c ($ x ) = 2e/hEjcos(2TT$ x /<S> ) of the SQUID, and E c , 
^bTjj -C Ej, we can neglect the driving term and ex- 
pand the term cos(ip) »1- ip 2 /2. As a result, we get a 
harmonic oscillator Hamiltonian with the mass given by 
the capacitance C, and the spring constant given by the 
Josephson inductance Lj — h/(2e)Ic- Now connecting 
the current fluctuations to the flux $ x = M(I + SI(t)), 
we can vary the inductance term, i.e., couple to x 2 of the 
harmonic oscillator. The coupling constant is given by 

Me 2 u> p . f n MI\ 

Q 1 1 = sin 2n 

yjJ 2h \ $ J ' 

where ui p — y/SEjEc/h is the plasma frequency of the 
SQUID at $ a = 0. With M = 20 pH, lo p =200 GHz, 
we get g — 10~ 3 . Hence, the current biased SQUID can 
be used for the detection of the third cumulant in the 
harmonic mode, but the parameters M and lo p need to 
be optimized to quite high values in order to obtain large 
enough coupling strength. 

Another way to use the SQUID in the harmonic mode 
would be to place it in a resonant circuit, and use the 
modulation of the Josephson inductance for the fluctu- 
ation measurement. In this case, the scheme would be 
similar to that presented in Subs. IVIDI but the coupling- 
would be to x 2 rather than p 2 . The coupling constant 
would be again given by Eq. IjVI E|) . but now u> p should 
be replaced with the resonance frequency of the circuit. 



measuring the third cumulant of current fluctuations in a 
short contact. When coupled to the driven non-Gaussian 
fluctuations, the static density matrix of these systems re- 
veals information on the frequency dependence of these 
fluctuations. The third cumulant can be read from the 
change in the effective temperature of these systems upon 
reversing the polarity of the bias across the contact. As 
the measured signal is inherently quantum, the order- 
ing of the current operators turns out to be important. 
Depending on the relative ratio of the frequency scales 
of the system, given by the voltage, eV/h, and the cir- 
cuit, ui c , the measured Fano factor for the third cumulant 
is either £„T„(1 - T„)(l - 2T n )/ Y, n T n (eV » hu; c , 
Eqs. m and |2T1» or 2Y /n T 2 {l - T„)/£„T„ {eV « 
hoj c , Eqs. (jTp|) and (J2^| ) ■ This slightly resembles the 
reasoning in Ref. 1461 where the measured Fano factor 
depends on "how far" the detector is placed from the 
current path. 
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VII. DISCUSSION 

We suggest to use the excitation and relaxation in 
quantum two-level systems or harmonic oscillators for 



The third-order 

h/4 
h/3 



Quite generally, the frequency dependent third cumu- 
lant considered in this manuscript satisfies 8 3 I(lui, U2) — 
5 3 I(—oji — u>2,u>2)- This is valid for a point contact at 
arbitrary temperature (Eqs. JSJ - 0) m the presence 
(5 3 I C from Eq. JSJ) or absence of the external circuit. 
Using this symmetry, it is straigthforward to show that 
h = — 14 and I3 = 

The integrals can be calculated at T — by 
evaluating (5 3 I(uj 7 ±cj ) — 5 3 I(~ w, ±u)q))/ui by parts, 
w € [0, min(w , \eV\/h — ujq)], oj € [min(wo, \eV\/h — 
uj ), max(w , \eV\/h — w )], w S [max(wo7 |eU|/?i — 
ujq), \eV\/h) and finally uo > \eV\/h. In the presence 



APPENDIX A: INTEGRALS REQUIRED FOR 
THE THIRD-ORDER RATES 



of the external circuit, described with Eq. ©, the re- 
sulting analytic expressions for are very long, even at 
zero temperature. However, they have rather simple lim- 
its depending on the relation between \eV\ and hujo, i.e., 
which of these scales gives the cutoff for the integrals. 
For all Hide, the integrals are finite only for \eV\ > O. 
Moreover, in both limits (but not generally), 1\ = 1-x- 



contributions to the rates depend on the integrals of the form 

6 3 I(lu,±lu ) - S 3 I{~LU,±LO ) 



Re 
Re 





= Re 


i 










= Re 


i 









du> 



dui 
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For \eV\ < hto c we get for T = 



Ji = eG(F 2 - F 3 )sgn(y) 



In 



|eV~| — Tiluq 



lant, J 3 /^,^) = G(u>i)G(uj2)G(— u)\ — u^i^eGV as 
was done in Rcf. 4. This yields 



|eV|ln 



\eV\ 



\eV\ — Ku)o 



h = —eGF 3 



2ttuj^luq 



(A2) 



(Al) 



For \eV\ ^> hui c , the integral can be evaluated by as- 
suming a frequency independent intrinsic third cumu- 



This result can be found fairly straigthforwardly via the 
residue theorem. 
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